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ABSTRACT 

We have made extensive observations of 35 distant slow (non-recycled) pulsars discovered in the 
ongoing Arecibo PALFA pulsar survey. Timing observations of these pulsars over several years at 
Arecibo Observatory and Jodrell Bank Observatory have yielded high-precision positions and mea- 
surements of rotation properties. Despite being a relatively distant population, these pulsars have 
properties that mirror those of the previously known pulsar population. Many of the sources exhibit 
timing noise, and one underwent a small glitch. We have used multifrequency data to measure the 
interstellar scattering properties of these pulsars. We find scattering to be higher than predicted along 
some lines of sight, particularly in the Cygnus region. Lastly, we present XMM-Newton and Chandra 
observations of the youngest and most energetic of the pulsars, J1856-i-0245, which has previously 
been associated with the GeV-TeV pulsar wind nebula HESS J1857+026. 

Subject headings: pulsars: general, pulsars: individual: PSR J1856+0245, surveys, ISM: structure; 
scattering 

1. INTRODUCTION 
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More than 200 pulsars are now known (Hobbs & 
Manchester 2013 ) . Most of these were discovered in bfind 
surveys by radio telescopes. In these surveys, short in- 
dividual telescope pointings, of order a few minutes, are 
made over grids of thousands of sky positions. In each 
pointing of a typical survey, the radio frequency signal 
is divided into hundreds of narrow spectral channels, 
each of which is square-law detected and sampled at a 
rate of order ten kilohertz. The resulting data blocks 
are searched for dispersed signals exhibiting either peri- 
odic or transient behavior expected of pulsars. For any 
newly discovered pulsar, survey data immediately yield 
approximate estimates of the pulsar's position (with the 
precision of the telescope primary beam size, typically 
between a few arc minutes and a degree) , along with the 
pulsar's period and dispersion measure. 

Follow-up timing observations of newly discovered pul- 
sars are needed to characterize more precisely the pulsar 
and its environment. Observations spaced over several 
months or longer decouple the influences of pulsar rota- 
tion and time-of-flight across the solar system, allowing 
measurement of pulsar spin-down rates and positions, the 
latter with sub-arcsecond precision. Such timing obser- 
vations can also detect orbital motion of pulsars in binary 
systems. Measured pulsar spin-down rates can be used 
to estimate the pulsar's age, magnetic field strength, and 
energy loss rate. Measured pulsar positions can be used 
to search for counterparts to the pulsar in other spec- 
tral bands, and phase coherent timing solutions can aid 
searches for pulsed radiation there. 

The Pulsar Arecibo L-Band Feed Array (PALFA) 
project is an ongoing deep pulsar survey of low Galactic 
latitudes being undertaken at the 305-m William E. Gor- 
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don Telescope at the Arec ibo Observatory ( Cordes et al 



2006 Lazarus et al. 2012). Among the survey goals are 



discovering millisecond pulsars suitable for use in pulsar 
timing arrays that will enable the detection of gravita- 
tional waves, discovering binary pulsars useful for tests of 
strong- field gravitation and determination of the neutron 
star equation of state, and discovering distant pulsars to 
better characterize the pulsar population and the inter- 
stellar medium on Galactic scales. 

In this paper, we describe several years of follow-up 
timing observations of 35 slow (non- recycled) pulsars dis- 
covered in the PALFA survey between 2004 and 2008. 
These are the first published observations of most of 
these sources, although a few of the discoveries have 
been announced p reviously eit her as conve ntional pul- 



sars ( [Cordes et ah [2006; Hessels et a l. [[ 2008[ ) or rotating 
radio transients (Deneva et al. 2009[ ). i'he timing ob- 
servations were rnade at both tne Arecibo Observatory 
and the Jodrell Bank Observatory. We give a synopsis of 
the PALFA survey in fJ2| We describe the observations 
in Sj3| We present timing models of the pulsars in ^ 
We give flux densities and spectral indices of the pulsars 
measured from the timing data in S|5[ We describe the 
pulse profiles and use them to characterize the scatter- 
ing of the pulsar signals by the interstellar medium in SM 
We describe dedicated X-ray observations of the young 
pulsar J1856-)-0245 in fJ7| We discuss our results in ^ 

2. THE PALFA SURVEY 
The paramete r s of the PALFA survey a r e disc ussed 



Cordes et al. 



(2006), van Leeuwen et al. (2006), and 



Lazarus et al. ([2012|). A detailed description of the early 



surve y observations will shortly be published (Swiggum 
2013 in preparation). We give a brief summary here. 

The survey area covers those portions of low Galactic 
latitude sky, |6| < 5°, visible with the Arecibo telescope. 
This includes Galactic longitude ranges 32° '^ i '^ 77° 
and 168° < ^ < 214°. Survey observations used the 
ALFA receiver, which provides seven simultaneous in- 
dependent beams on the sky. Data were collected over 
100 MHz passbands centered on 1420 MHz in each re- 
ceiver beam. Each passband was processed by a Wide- 
band Arecibo Pulsar Processor 3-level autocorrelator sys- 
tem (see also §[3]). Autocorrelation spectra with 256 lags 
were recorded at intervals of 64 /is, with orthogonal po- 
larizations summed before recording. The typical dwell 
time for a telescope pointing was 268 s. (More recent 
PALFA survey observations, made after the discoveries 
discussed in this paper, have used upgraded spectrome- 
ters, providing 960 channels across 320 MHz passbands.) 

Previous surveys of the low Galactic latitude sky vis 



ible from Arecib o (Hulse & Taylor 1975 Stokes et al 



[1986! Nice et aL[ 19SS ) used lower radio frequencies and 
were constrained by the narrow bandwidths available on 



were constrainecl by the narrow bandwidths available on 
the telescope before its upgrade in the mid-1990s. In 
contrast, the PALFA survey's relatively high radio fre- 
quency of 1420 MHz, relatively large bandwidth, and fine 
spectral resolution give it sensitivity to pulsars over a far 
larger volume of the Galaxy than ever before, particu- 
larly to those pulsar s with short spin peri ods and high 
dispersion measures ( Crawford et al.[p012 ). 

The bulk of the pulsars observed lor this paper were 
discovered by the "quick look" pipeline, a low-resolution 
periodicity search, which ran in near-real-time as data 



were collected ( Cordes et al.[[2006 1 . The remainder were 
discovered either by tuU-resolution periodici ty searches or 



by se arches for dispersed transient signals ( [Cordes et al. 
2006 ; [Deneva et aL[[2009| ) . 



While the present paper focuses on isolated long-period 
pulsars, the PALFA survey has also found numerous bi 
nary and millisecond pulsars (e.g., [L orimcr ct al. 2006 



Champion etaL[[2008l IKnispel et ^ ^010, |2011, |Deneva| 



I 



et al.[[2012| [Crawford et al.[[2012| [Allen et al.[[2013[ ). 

3. TIMING OBSERVATIONS 

We observed all thirty-five pulsars in this paper at the 
Arecibo Observatory. In addition, we observed seventeen 
of the stronger sources at Jodrell Bank Observatory. De- 
tails of observations for individual sources (observatory 
used, cadence, time span of observation) can be inferred 
from plots presented in Sj4j below. Although several of 
the pulsars in this paper were discovered through detec- 
tion of individual pulses, all timing for this paper was 
done by averaging pulsar time series over many pulse 
periods. 

Observations at Arecibo Observatory used the "L- 
wide" receiver, with system equivalent flux density of 
approximately 2.4 to 3.0 Jy, depending on zenith angle, 
frequency, and sky temperature. Data were collected us- 
ing four Wideband Arecibo Pulsar Pr ocessors (WAPPs ) 
3-level autocorrelation spectrometers ( Dowd et al.|2000" ). 
The spectrometers operated in parallel to cover tour re- 
ceiver subbands, centered at frequencies of 1170, 1370, 
1470, and 1570 MHz, and with 100 MHz bandwidth in 
each subband. Within each subband, 256-lag autocor- 
relations were accumulated and recorded at intervals of 
256 fis. Self- and cross-products of two circular orthog- 
onal polarizations were recorded, although only the self- 
products were used for the present work. A typical ob- 
servation duration was 2 to 3 minutes. A pulsed radio 
frequency reference signal of known strength was injected 
into the receiver and recorded prior to the observation of 
each pulsar in order to calibrate the receiver ampliflca- 
tion, as needed to balance the two orthogonal polariza- 
tions and to measure absolute fiux densities. Absolute 
flux calibration used measurements of the reference sig- 
nal strength provided by the Observatory. 

The recorded spectral time series were dedispersed and 
folded off-line u sing a modifled version of the SiGPROC 
software packag q^ Data from each subband were pro- 
cessed separatelyTObservations were processed in subin- 
tervals of 30 seconds, so a typical 2-3 minute observa- 
tion in four subbands yielded 16 to 24 folded pulse pro- 
flies. Each of these profiles was cross-correlated with a 
reference profile to produce a time of arrival (TOA) us - 



ing standard techniques (e.g., Lorimer fc Kramer[[2005 ) 
Thus the Arecibo data set for any given observation of 
a source consisted of order two dozen TOAs. The same 
reference profile was used to process data from all four 
subbands, even for pulsars with significant profile evo- 
lution. The TOA sets for some weak pulsars included 
a significant number of TOAs which did not fit timing 
models (described below), because of signal-to-noise ra- 
tios too low to robustly detect the pulsar signal. The 
data sets were edited using manual and semi-automated 
methods to remove obviously spurious TOAs. 



■^^ http://sigproc.sourceforge.net 
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Observations at Arecibo Observatory were made at 
roughly six- week intervals between 2006 March and 2010 
March except for a gap between 2007 March and 2007 
December due to extended telescope maintenance. A 
few additional observations were made in 2012 Septem- 
ber. Within an observing session, any given pulsar was 
observed either on a single day or a on pair of adjacent 
days. Occasionally pulsars were observed multiple times 
within a day to resolve pulse numbering ambiguities. 

Observations at Jodrell Bank used a dual-polarization 
cryogenic receiver on the 76-m Lovell telescope, with a 
system equivalent flux density of 25 Jy on cold sky. Data 
were processed by a digital filterbank which covered the 
frequency band between 1350 MHz and 1700 MHz with 
channels of 0.5 MHz bandwidth. Observations were typ- 
ically made with a total duration of between 10 minutes 
and 40 minutes, depending upon the discovery signal-to- 
noise ratio. Data were folded at the nominal topocen- 
tric period of the pulsar for subintegration times of 10 
seconds. After inspection and "cleaning" of any radio- 
frequency interference, the profiles were dedispersed at 
the nominal value of the pulsar dispersion measure. Ini- 
tial pulsar parameters were established by conducting lo- 
cal searches in period and dispersion measure about the 
nominal discovery values and finally summed over fre- 
quency and time to produce integrated profiles. TOAs 
were obtained after matching with a standard template 
and processed using standard analysis techniques with 
PSRTiMEFland TempcEB 

Observations at Jodrell Bank were made between 2005 
June and 2012 April. Observations of a given source were 
made at intervals of 1 to 2 weeks. The observation of any 
given source yielded a single TOA on a given day. 

4. TIMING 

The TOAs from Arecibo and Jodrell Bank (where ap- 
plicable) were combined into a single data set for each 
pulsar, which was fit to a standard timing model using 
X^ minimization techniques implemented in the Tempo 
software package. The basic timing model for each source 
was parameterized by right ascension, a; declination, S; 
pulse phase; rotation period, P; rotation period time 
derivative, P; and dispersion measure, DM. Additional 
parameters were fit to model timing noise in some pul- 
sars, and a glitch in one pulsar, as described below. The 
best-fit parameters for each pulsar are listed in Table [T] 
Parameter uncertainties listed in the table are twice the 
formal uncertainties calculated in the timing model fit. 
Spin periods and period derivatives are plotted in Fig- 
ure[T] Values of DM given in this table do not includ e cor - 
rections for scattering or intrinsic profile evolution (S 6.2 ). 

For pulsars observed at Jodrell Bank, an arbitrary time 
offset was allowed between the Arecibo and Jodrell Bank 
TOAs. The timin g analysis used the JPL DE405 solar 
system ephemeris ( |Standish|1998l ) and the TT(BIPMll) 
time scaler^ Because of the prevalence of timing noise 
(see belowJT unweighted fits were used. 

Residual pulse arrival times from fits of each pul- 
sar to the basic timing model are shown in Figures |2|- 



[6j Timing noise is evident in many sources. To 
measure accurate pulsar parameters, TOAs for these 
sources were "whitened" by including a Taylor expan- 
sions of the pulsar period in the timing model, P(t) — 
E"=o(l/^0-P^"Hi-*o)", where t is time, to is the period 
epoch, and P'"' = [d^P/dt'^]tg. For each pulsar, mod- 
els with rifit from 1 (the basic timing model) through 10 
were tried, and an appropriate value of nm was chosen 
based on the standard deviations of the residuals and 
on visual inspection of the residuals to insure there were 
not groups of correlated nonzero residuals. The chosen 
values of ngt for each pulsar are listed in Table [T] 

For four pulsars, designated as "N" in the ngt column 
of Table [T] polynomial models with ngt as high as 10 did 
not yiela white residuals. An alternative method was 
used to estimate the value of each timing parameter for 
these pulsars. For each pulsar, the timing model (includ- 
ing tenth order polynomial in P(i)) which minimized x^ 
was found. Each parameter was then individually per- 
turbed (allowing other parameters to freely vary in the 
fit) to find the parameter values which resulted in a three- 
fold increase in the standard deviation of the residuals. 
The resulting range of parameter values was used to set 
the uncertainties presented in Table [l] We view this as 
a highly conservative approach to estimating the uncer- 
tainties. 

Because many of the pulse periods were analyzed as 
high-order Taylor series around the period epochs, the 
values of P and P listed in Table [ij are the best estimates 
of the values at the period epoch; they are not averages 
over the length of the data set. 

In Table [2] we present several quantities derived from 
the timing model parameters. Estimates of spin-down 
age, tg, surface magnetic field strength, B, and spin-down 
energy loss rates, E, were made using conventional for- 
mulas (Lorimer & Kramer 2005) and are given in loga- 
rithmic iornrTn~tIie~ta5Ie7'^Positions in Galactic coordi- 
nates are based on the positions in equatorial coordinates 
in Table [3] Distances are estimated from dispersion mea- 
sured in Table [1] using the NE2001 mode l of the distribu- 
tion o f ionized material in the Galaxy ( jCordes fc Lazio 
2002|). The positions of the pulsars projected onto the 
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Galactic plane are shown in Figure [7] 

We detected one glitch, in PSR J1947-f 1957, as evinced 
by the peak in the residual pulse arrival times for 
this pulsar in Figure |6] The glitch event was on 
MJD 55085.5 ± 2.6 and was well fit by an instanta- 
neous, permanent change in rotation frequency of Az/ — 
(1.02 ± 0.03) X 10^^ s^^, which is a fractional change of 
Aiy/iy = (1.61 ± 0.04) x 10"^. This is a typ ical value 
for a small pulsar glitch (e.g., Espinoza et al. 2011). 
Other pulsars in our data set may have had glitch-like 
events; in particular, see the residual arrival times of 
PSRs J1905+0902 and J1916-H12245 in Figure [3] How- 
ever, due to the small sizes of the events and the coarse 
sampling of the time series, these events cannot easily be 
distinguished from smoothly-varying timing noise. 

5. FLUXES AND SPECTRAL INDICES 

Flux densities and spectral indices of the pulsars are 
given in Table |4] These were measured from the cali- 
brated pulse profiles from the Arecibo observations by 
the following procedure. (For logistical reasons, only 
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data from some observing dates were used.) Each pro- 
file was generated from a 30-second observation in one of 
the four 100-MHz subbands, as described above. The 
flux density was measured in each profile. All flux 
density measurements for a given subband in a given 
day were combined; the resulting daily average values 
were then combined to find a global average for that 
subband. Finally, the four average flux densities from 
the four subbands were fit to a function of the form 
S{f) = S'i4oo(/mhz/1400)", with flux density 51400 and 
spectral index a as free parameters in the fit. The 
best-fit parameters for each pulsar are given in Table |4l 
along with the reduced chi-squared values, x^, of the 
fits. Although the formal uncertainties from the fits are 
reported in the table, they should be interpreted with 
caution, both because they arise from fits with only two 
degrees of freedom (four data points and two fit parame- 
ters) and because scintillation could cause apparent, non- 
intrinsic variations in pulsar flux densities from epoch 
to epoch. Table [4] also includes luminosity estimates, 
Liioo — Si4QocP, using the distances derived from DMs 
as reported in Table [3J 

6. PROFILE MORPHOLOGIES AND SCATTERING 
PROPERTIES 

6.1. Profiles 

Profiles of each pulsar in each of the four observing 
subbands are given in Figures 8pT These average pro- 
files were generated by coherently adding all Arecibo 
data profiles for a given pulsar that yielded good TOAs. 
The alignment of the profiles was made using the timing 
solutions described in fj4] By and large, these profiles 
are typical of those seen in the previously known pul- 
sar population. Some evidence of profile evolution over 
frequency is seen (e.g., relative component strengths of 
PSR J1948-f2551). 

One source, J1909-I-0749, has a strong interpulse. Of 
its two pulse components, we have somewhat arbitrarily 
designated the shorter, broader component (near the cen- 
ter of its profiles in Figure pi) as the main pulse and the 
taller, narrower component (near the end of its profiles 
in Figure ^ as the interpulse. The interpulse trails the 
main pulseby ^ 0.484 of a pulse period (174°). The ratio 
of integrated flux density for the main pulse over that of 
the interpulse in each of the Arecibo observing subbands 
is 1.28, 1.21, 1.02, and 1.36, at 1170, 1370, 1470, and 
1570 MHz, respectively. We do not know of a physical 
explanation for the dip in this ratio at 1470 MHz. 

The profiles in Figures [8] and |9] were used to compute 
pulse widths at half-maximum, W50, for each pulsar in 
each observing band. The values are given in Table [5j 
both in milliseconds and as a fraction of a pulse period? 

6.2. Scattering 

The profiles of several pulsars show evidence of 
frequency-dependent pulse broadening, indicative of mul- 
tipath scattering in the interstellar medium. The ob- 
served profiles result from the convolution of pulse- 
broadening functions (PBFs) with the intrinsic pulse 
shapes. To quantify the scattering for each pulsar's sig- 
nal, we used least-squares techniques to simultaneously 
fit PBFs to the profiles of the pulsar in each of the four 
observing subbands. While PBFs can take on a range 



of shapes dependent on the distribution of scattering 
material along the line of sight and on the wavenumber 
spectrum of irregularities, for simplicity we adopt a one- 
sided exponential form with a characteristic e~^ broad- 
ening time that scales as Tissif) = '''ighz(//1 GHz)^**. 
We modeled the intrinsic profile of each pulsar as a sum 
of between one and five Gaussian components. The posi- 
tion and width of each component was held fixed across 
all four subbands, but the amplitudes were free to vary. 
The fit included scattering time, tighz! an arbitrary off- 
set in DM, ADM, to account for biases in the DM used 
to generated the profiles (since the dedispersion used for 
those profiles did not account for scattering); the posi- 
tion and width of each component; and the amplitude of 
each component in each subband. Thus a total of 8 to 
32 parameters were free to vary in the fit for each pulsar, 
depending on the number of Gaussian components in the 
model intrinsic profile. Some profiles had noisy baselines, 
presumably due to radio frequency interference; a poly- 
nomial (up to twentieth order) was included in the model 
fit to such profiles. 

This model proved satisfactory for all pulsars in our 
sample. Examples of fitted pulse profiles are in Figure [12} 
The results of all the fits are given in Table |3j which lists 
TiGHz; ADM; the number of components in the profile 
model, ncomp; and the goodness-of-fit statistic, xl- Scat- 
tering was detectable in a large number of the sources, 
and the time scales were determined with precision of 
a few percent for the most highly scattered sources. 
The fits are characterized by reduced chi-squared values 
around xt = 1.0±0.5 for most sources. Even though this 
method fits the pulse profile well, there is an uncertainty 
in assuming Gaussian-shaped components in the intrinsic 
profile models. For a multi-component model combined 
with a small amount of scattering, the second or third 
component is covariant with scattering time. We are not 
able to resolve the ambiguity with the current data. Ob- 
servations at lower frequency (and hence higher scatter- 
ing) and measurements of scintillation bandwidths could 
help resolve these covariances. 

We have not tested for departures from the assumed 
frequency dependence tjss oc /~^ or for the dependence 
of the results on alternative forms for the PBFs, such as 
those that result from electron-density variations having 
a Kolmogorov wavenumber spectrum and distributed in 
a thin sc reen or an extended medium (Lambert & Rick- 
1999|). Anisotropic scattering also presents alterna- 



ett 



tive PBF forms as does any variation of the scatte ring 
transverse to the line of sight ( Cordes fc Lazio||2001[ ). 

A comparison between our measured scattering time 
scales for pulsars in the inner Galaxy (30° < I < 80°) 
and the predictions of t he NE2001 electron-density model 
(Cordes & Lazio 2002) is given in Figure 13 The pulse 
broadening times m our sample of pulsars displays the 
same amount of variation at a given value of DM as is 
seen in oth er objects (e.g., 'Cordes & Lazio 2002 Bhat 



et al.||2004[ ). While the NE2001 electron-density moHer 
generally underpredicts the pulse broadening times for 
the pulsars in our sample with the largest pulse broad- 
ening, the broadening times are consistent with the ob- 
served ranges from the larger sampl e of objects used t o 



construct NE2001 and analyzed in Bhat et al. (2004) 



The underprediction is especially notable for objects m 
the Cygnus direction (as seen at the higher Galactic Ion- 
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gitudes in Figure 13 ). This is to be expected because the 
Cygnus superbubole comprises discrete HII regions tha t 



are known to produce excess scattering ( Fey et aTjigSQ ) 
All but one of the pulsars (PSR J2UU5+3547J have 
DMs that can be accommodated by the NE2001 model, 
but some of the DM values and pulse broadening times 
are undoubtedly enhanced by intersection of the line of 
sight with particular electron-density enhancements not 
included in NE2001. A new electron density model is 
now under development that will include the DM values 
and scattering measurements discussed here. 

The scattering measurements imply that searches for 
millisecond pulsars along low-Galactic-latitude lines of 
sight in our observing region of 32° < i < 11° will 
be limited by scattering at very large distances. Pre- 
suming a scaling law of riss oc /~^, the median mea- 
sured value TiGHz ~ 4 ms among our pulsars scales to 
■'"1.42GHz ~ 1 ms. Scattering of this magnitude signifi- 
cantly curtails sensitivity to the fastest-spinning pulsars, 
especially those with periods of order 1 millisecond. (It 
is important to note that there is still a large volume of 
space, likely extending to distances of several kiloparsecs, 
for which scattering is not significant, and within which 
there is promise for a large number of millisecond pulsars 
to be detected by surveys such as PALFA.) 

7. X-RAY OBSERVATIONS OF PSR J1856-h0245 

The 81-ms pulsar J1856-I-0245 has the smallest char- 
acteristic age, 21 kyr, and largest spin-down luminosity, 
E = 4. 6 X 10^^ergs~\ of the pulsars observed for this 
paper. Hessels et al. (2008) found this pulsar to be coin- 
cident with the extended ieV source HESS J1857-I-026 
and noted faint X-ray emission in the pulsar's vicinity 
in archival ASCA data. They concluded that HESS 
J1857-I-026 was plausibly a pulsa r wind nebula (PWN) 
associated with PSR J1856-t-0245. [Rousseau eral| ( |2012 1 
examined 36 months of Fermi Large Area i'elescope 7-ray 
data toward this pulsar. They detected emission coinci- 
dent with the pulsar and extended TeV source, but they 
did not detect pulsations. They interpreted the unpulsed 
emission from the region as a putative GeV PWN, and 
performed multi-wavelength (keV-TeV) modeling of the 
spectral energy distribution. Here we describe dedicated 
XMM-Newton and Chandra X-ray observations of this 
region, which were obtained to verify the PWN interpre- 
tation. These observations were used in lRousseau et al.1 
(2012) to place a limit on the flux of any extended X-ray 
emission. Here we present for the first time the X-ray 
spectrum of the pulsar and further discuss its associa- 
tion with the GeV-TcV emission. 



7.1. XMM-Newton 

XMM-Newton observed the field of PSR J1856-fG245 
on 2008 March 27-28 (ObsID 0505920101) for a total 
of 55.3 ks. Data were acquired from only the European 
Photon Imaging Camera (EPIC) pn detector, operating 
in full-frame imaging mode with the thick optical filter 
in place. Unfortunately, due to tests being conducted si- 
multaneously on the MOSl and M0S2 detectors, no data 
are available from those instruments. The EPIC pn field 
of view was centered close to the timing position of PSR 
J1856-I-0245 and also included the centroid of the TeV 



pn data were reprocessed using the Epchain pipeline 
in SASJ^The data were subsequently screened for in- 
stances of strong soft-proton fiaring. The filtering pro- 
cedure left 30.2 ks of useful exposure time. The Refiec- 
tion Grating Spectrometer (RGS) data provided no use- 
ful source information and are not included in our anal- 
ysis. The SAS source detection task edetect.chain 
found a source coincident with the radio timing position 
of PSR J1856-h0245, XMMU J185650.8-h024545. 

Photons were extracted from a circular region of ra- 
dius 10" centered on the pulsar, which encloses ~52% of 
the total energy. This relatively small extraction radius 
was chosen to minimize contamination from the nearby 
(~ 22") source XMMU J185651. 8-1-024530. An estimate 
of the background was obtained from several source-free 
regions around the pulsar position. The spectral analy- 
sis was performed using XSPECJ^ 12.5.0.5. The source 
spectrum was grouped to ensure at least 15 photons per 
energy bin. 

The X-ray counterpart of PSR J1856-h0245 shows a 
hard spectrum, most likely because it is distant and heav- 
ily absorbed. We fit the spectrum of the source with 
both absorbed power-law and thermal models. Given 
the low photon count, both models describe the data 
adequately well, though an absorbed blackbody model 
gives a temperature of kT w 1.5 — 3keV, too hot to be 
the surface temperature of the neutron star. For this 
reason, the power law model is preferred. For young 
pulsars, such emission is typically due to particle accel- 
eration in the pulsar magnetosphere. The best-fit pa- 
rameters for this model are only weakly constrained due 
to the limited number of photons with absorption TVjj 
(4.9 ' ' 



-1-3. 2^ 
-2.4^ 



X 10 cm "^ and photon index F = 1.2; 



-14 



-0.8 rp, 

-0.4- -'^"'= 



nebula HESS J1857+026 (Figure 14). The raw (ODF) 



implied unabsorbed flux is ~8.3 x 10 "" erg cm 
(2-10 keV), which for a distance of 9 kpc (Table [3]) cor- 
responds to a luminosity of Lx ^ 8 x 10'^^ erg s . 

The 73.4 ms time resolution of the EPIC pn full-frame 
exposure and the 3.2 s frame time of the ACIS-I observa- 
tion preclude the search for X-ray pulsations from PSR 
J1856-f0245. 

7.2. Chandra 

PSR J1856+0245 was observed with Chandra ACIS-I 
for 39 ks on 2011 February 28 (ObsID 12557). These 
data were recorded in VFAINT and Timed Exposure 
(TE) modes and were analyzed using CIAGrj version 
4.3.1 with CALDB 4.4.3. PSR J1856-f0245 is clearly de- 
tected as a poin t sourc e (Figure 14). As also discussed in 
Rousseau et al. (2012), there is no evidence for extended 
emission surrounding the pulsar. Subtracting a simu- 
lated point spread function using CHARlj^and MARxj^ 
from the observed point source yields no excess counts 
that could arise due to a compact PWN. Extracting the 
emission from an annular region with inner and outer 
radii of 2" and 15", respectively, we find an upper limit 

^'' The XMM-Newton SAS is developed and maintained by the 
Science Operations Centre at the European Space Astronomy Cen- 
tre and the Survey Science Centre at the University of Leicester. 

^* http : //heasarc . nasa . gov/docs/xanadu/xspec/index . html 

^^ Chandra interactive Analysis of (J bservations (l^'ruscionc et al.| 
'2006), http : //cxc . harvard . edu/ciao/ 

•^" I'he Chandra Kay Tracer, http : //cxc ■ harvard. edu/soft/| 
|ChaRT/cgi-bin/www-saosac. cgi 

'^^ fhttp: //space. mit .edu/cxc/marx/ index. html 
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on the unabsorbed flux of 5 x 10^'' erg s~^ cm^^ (1-10 
keV, 3(7 confidence), corresponding to a luminosity of 
<5 X 10^^ erg s~^. In deriving this limit we assumed a 
power-law index of F — 1.5 and an equivalent neutral 
hydrogen column depth iVn = 4 x 10^^ cm~^. The latter 
is somewhat higher than would be predicted by a sim- 
ple em pirical correlation between A^h and DM (He et al. 
2013|, so can be regarded as conservative. Note that the 



outer extraction radius was chosen based on the X-ray 
PWN for pulsars with comp arable E by scaling their an - 
gular size with the distance ( Kargaltsev fc Pavlov|20d8 ), 
assuming a distance of 9 kpc ('i'abie|3|. 



7.3. Association between PSR J1856+0245 and HESS 
J1857+026 



Hessels et al. (2008)) showed that PSR J1856-I-0245 
is coincident witn the faint ASCA X-ray source AX 
J185651-I-0245 and suggested that the latter might rep- 
resent an X-ray PWN associated with the pulsar. Our 
XMM-Newton and Chandra observations confirm X-ray 
emission from the pulsar itself, but find no evidence of 
surrounding diffuse emission. While the discovery of an 
X-ray PWN — especially one which shows a morphology 
indicative of an association with the TeV PWN (e.g., 

would strengthen the link be- 



Aharonian et al. 2005 



tween PSR J 1856-^0245 and HESS J1857-^026, the lack 
of an observable X-ray PWN does not eliminate the pos- 
sibility of association. Given the distance of 9 kpc in- 
ferred from the pulsar's DM, it is possible that the X-ray 
PWN is simpl y too faint to be det ected in the current 
data. Indeed, Rousseau et al. (2012) show that it is pos- 
sible to derive a self-consistent model of the spectral en- 
ergy distribution that accounts for the low X-ray nebula 
flux. 

We have taken a critical look at other X-ray sources de- 
tected in our XMM-Newton and Chandra observations, 
and find no evidence for extended sources of a possible 
PWN nature (Figure 14). These observations do not 
cover the full extent of the TeV nebula, but they do 
include the source centroid. Furthermore, the PALFA 
survey has covered the region of the HESS TeV emission 
and as such any additional, undiscovered, young pulsar in 
the region would have to be quite faint, 5*1400 ^ 0.2 mJy 
(depending on the duty cycle, dispersion measure, and 
scattering). We thus conclude that PSR J1856-I-0245 re- 
mains the most viable association for HESS J1857-I-026. 

8. DISCUSSION 

We have presented timing properties and pulsar pro- 
file morphologies of thirty-five pulsars discovered in the 
PALFA survey. Due to this survey's depth and its sen- 
sitivity to high-DM pulsars, these pulsars are relatively 
distant, with a median distance of 7.1 kpc compared to 
the median distance of 4.2 kpc of Galactic pulsars listed 
in the ATNF pulsar catalog (Figure [TJ Table [3|. 

Despite their large distances, the properties of these 
pulsars are similar to those of the previously known 
non-recycled Galactic pulsar population. We used 
Kolmogorov-Smirnov tests to compare the distributions 
of spin-down ages, magnetic field strengths, and spin- 
down energies of the pulsars in this paper with those 
of the known pulsar population. We excluded recy- 
cled pulsars from this comparison by imposing a cut- 



off P > 10^"'^^. We found no significant differences be- 
tween the two populations. (A more thorough analysis 
of the population of pulsars discovered by PALFA, using 
more powerful statistical techniques and considering the 
details of the telescope search pointings, is beyond the 
scope of the present paper.) 

Four of the sources, PSRs J1856-h0245, J1909-h0749, 
J1928-hl746 and J1934-h2352, have spin-down ages t^ < 
10^ yr. Pulsars of such young age are often found in 
supernova re mnants. Howe ver, Green's Supernova Rem- 
nant catalog (Green 2009) lists no remnants coincident 
with the positions of these pulsars. The remnant catalog 
is known to be incomplete, however, and dedicated radio 
interferometric observations of these objects in particular 
to look for remnants could be fruitful. 

Young and energetic radio pulsars often have coun- 
terparts at other wavelengths, such as the X-ray source 
at the position of PSR J1856-I-0245 described above. 
We have checked the ROSAT all-sky survey for X-ray 
point sources at the positions of all the pulsars having 
E > 2 X 10'^^ erg s~^, but we have found no counter- 
parts. 

As discussed in 34) many of the pulsars in thispa- 



per show deviations from steady spin-down (Figures 2 



6) 



The spin evolution of a pulsar is conventionally quanti- 
fied by a braking index, n, defined through v (x v^ , where 
V = 1/P is the spin frequency and u its time deriva- 
tive. For a pulsar undergoing energy loss due to magnetic 
dipole rotation, n = 3. In principle, the braking index for 
a pulsar can be calculated from n = vi) /v^ if the second 



time derivativ e of spin frequency, i), is known (Lorimer & 
(Kramer 2005 ) . In practice, such meas urements has been 
achie ved tor only a few pulsars (e.g., iLivingstone et al. 
2005[ and references therein). Most measurements of v 
and inferred values of n have been much larger in magni- 
tude than pr edicted by steady spin-down models. )Hobbs| 



I 



et al. (2010) analyzed long-term pulsar timing data of 
hundreds 01 pulsars and found that pulsars with spin- 
down ages less than 10^ yr all had positive values of n, 
whereas older pulsars had a mix of positive and negative 
values. They interpreted this as evidence that young pul- 
sars are recovering from glitches that occurred prior to 
the starts of the data sets, while older pulsars are domi- 
nated by stochastic processes. 

Significant values of i^, and hence n, can be mea- 
sured for about half of our sources (all those for which 
nfit > 1 in Table fTl). The four pulsars with spin- 
down ages less than 10*^ yr all show positive values of 
n: n ^ 24.5, 11.3, 22.9, and 12.3 for PSRs J1856-h0245, 
J1909-I-0749, J1928-H1746, and J 1934-h2352, respec - 
tively, consistent with the findings of'Hobbs et al. (2010). 
The older pulsars show a mix of positive and negative val- 
ues of n. All measured values of n are sufficiently large in 
magnitude to suggest that they are not associated with 
steady spin-down of these pulsars. 

The dispersion measurem ents and scattering measure- 
ments presented in f|4]and ^6.2 are a promising start of 
the detailed study ofthe ionized interstellar medium in 
this range of Galactic longitudes. The scattering mea- 
surements could be greatly enhanced by collecting and 
analyzing profiles of these sources at longer wavelengths. 

The 35 pulsars in this paper are just the tip of the ice- 
berg. As of this writing, the ongoing PALFA survey had 
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discovered more than 100 pulsars. Study of these objects 
will enhance our knowledge of the pulsar population and 
the interstellar medium in this sector of the Galaxy. 
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Table 1 

Timing Parameters of 35 Pulsars^ 



PSR 


Previous 


a 


5 


Epoch of P 


P 


P 


DM 


"fit 




Reference*^ 


(J2000) 


(J2000) 


(MJD) 


(s) 


(10-15) 


(pccm"'') 




J0540+3207 
J0628+0909 
J1848+0351 
J1850+0423 
J1855+0205 


1 
1,2 


05:40:37.116(1) 

06:28:36.183(5) 

18:48:42.24(1) 

18:50:23.43(1) 

18:55:42.046(5) 


32:07:37.3(1) 
09:09:13.9(3) 
03:51:35.7(3) 
04:23:09.2(4) 
02:05:36.4(2) 


54780.000 
54990.000 
54620.000 
54600.000 
54510.000 


0.5242708632083(5) 

1.241421391299(3) 

0.191442663494(2) 

0.290716217162(3) 

0.2468167699732(9) 


0.44811(2) 

0.5479(2) 

0.0673(1) 

0.0914(2) 

0.06466(7) 


61.97(4) 

88.3(2) 

336.6(4) 

265.8(4) 

867.3(2) 


1 
1 
1 
1 
1 


J1856+0245 
J1858+0346 
J1904+0738 
J1905+0902 
J1909+0641 


3 

1 
1 
2 


18:56:50.9(3) 

18:58:22.36(2) 

19:04:07.533(1) 

19:05:19.535(2) 

19:09:29.052(4) 


02:45:47(9) 

03:46:37.8(8) 

07:38:51.69(4) 

09:02:32.49(8) 

06:41:25.8(2) 


54930.000 
54510.000 
54760.000 
54570.000 
54870.000 


0.08090668906(4) 

0.256843797950(5) 

0.2089583321216(2) 

0.2182529126846(9) 

0.741761952452(6) 


62.117(4) 

2.0401(4) 

0.410898(7) 

3.49853(8) 

3.2239(7) 


623.5(2) 

386(1) 

278.32(8) 

433.4(1) 

36.7(2) 


N 
1 
1 
3 
1 


J1909+0749 
J1916+1225 
J1917+1737 
J1919+1314 
J1919+1745 


2 


19:09:08.2(2) 

19:16:20.045(1) 

19:17:23.976(5) 

19:19:32.986(5) 

19:19:43.342(4) 


07:49:32(5) 

12:25:53.94(4) 

17:37:31.6(1) 

13:14:37.3(1) 

17:45:03.79(8) 


54870.000 
54570.000 
55070.000 
54790.000 
55320.000 


0.23716129322(9) 

0.227387488792(2) 

0.334725226760(4) 

0.571399860595(2) 

2.081343459724(9) 


151.920(9) 

23.4516(2) 

0.324(1) 

3.78771(6) 

1.7050(4) 


539.36(5) 

265.31(3) 

208.0(2) 

613.4(2) 

142.3(2) 


N 
8 
1 
1 
1 


J1924+1631 
J1928+1746 
J1934+2352 
J1938+2010 
J1940+2337 


1 


19:24:54.88(3) 

19:28:42.553(1) 

19:34:46.19(2) 

19:38:08.34(2) 

19:40:35.486(7) 


16:31:48.8(6) 

17:46:29.62(3) 

23:52:55.9(3) 

20:10:51.7(4) 

23:37:46.5(2) 


54690.000 
54770.000 
54650.000 
54940.000 
54560.000 


2.9351864592(2) 

0.0687304058253(2) 

0.17843152366(2) 

0.68708185664(2) 

0.546824157964(7) 


364.212(4) 

13.190700(6) 

130.699(2) 

3.399(1) 

76.7965(2) 


518.5(9) 

176.68(5) 

355.5(2) 

327.7(8) 

252.1(3) 


2 
4 
8 
4 
2 


J1941+2525 
J1946+2535 
J1947+1957 
J1948+2551 
J1949+2306 




19:41:20.80(1) 

19:46:49.10(5) 

19:47:19.435(4) 

19:48:17.579(1) 

19:49:07.32(1) 


25:25:05.3(2) 

25:35:51.5(9) 

19:57:08.39(9) 

25:51:51.95(2) 

23:06:55.5(2) 


54920.000 
54810.000 
55190.000 
54930.000 
54910.000 


2.30615269700(2) 

0.51516705131(6) 

0.157508542541(4) 

0.1966268136138(4) 

1.31937338530(1) 


160.8349(5) 

5.641(5) 

0.5226(1) 

9.02244(2) 

0.1231(3) 


314.4(4) 

248.81(4) 

185.8(2) 

289.27(5) 

196.3(5) 


1 

N 
1 
3 

1 


J1953+2732 
J2005+3547 
J2006+3102 
J2007+3120 
J2009+3326 


1 


19:53:07.80(2) 

20:05:17.493(9) 

20:06:11.0(3) 

20:07:09.012(6) 

20:09:49.16(2) 


27:32:48.2(4) 

35:47:25.4(1) 

31:02:03(4) 

31:20:51.54(9) 

33:26:10.2(3) 


54570.000 
54540.000 
54800.000 
54920.000 
54430.000 


1.33396499366(2) 

0.615033894871(4) 

0.16369523645(9) 

0.60820531457(1) 

1.43836860189(2) 


1.773(1) 

0.2808(3) 

24.872(9) 

15.604(1) 

1.468(2) 


195.4(9) 

401.6(3) 

107.16(1) 

191.5(2) 

263.6(8) 


1 
1 
N 
9 
1 


J2010+2845 
J2010+3230 
J2011+3331 
J2013+3058 
J2018+3431 


1 
1 

1 


20:10:05.068(2) 

20:10:26.50(1) 

20:11:04.926(2) 

20:13:34.253(3) 

20:18:53.196(2) 


28:45:29.17(3) 

32:30:07.3(2) 

33:31:24.64(3) 

30:58:50.69(5) 

34:31:00.51(2) 


54810.000 
54430.000 
54490.000 
54600.000 
54760.000 


0.5653693451984(6) 

1.44244752060(1) 

0.931733093470(4) 

0.2760276934460(5) 

0.3876640866385(8) 


0.09072(2) 

3.616(1) 

1.7857(1) 

0.15194(4) 

1.83649(2) 


112.47(8) 

371.8(5) 

298.58(6) 

148.7(1) 

222.35(7) 


1 
1 
2 
1 
2 


^ Figures in par 
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Table 2 




Quantities Derived from Flotation Parameters 




of Table [T] 




PSR 


log ts 


logB 


\ogE 




(logyr) 


(logG) 


(logerg/s) 


J0540+3207 


7.27 


11.7 


32.1 


J0628+0909 


7.56 


11.9 


31.0 


J1848+0351 


7.65 


11.1 


32.6 


J1850+0423 


7.70 


11.2 


32.2 


J1855+0205 


7.78 


11.1 


32.2 


J1856+0245 


4.31 


12.4 


36.7 


J1858+0346 


6.30 


11.9 


33.7 


J1904+0738 


6.91 


11.5 


33.2 


J1905+0902 


5.99 


11.9 


34.1 


J1909+0641 


6.56 


12.2 


32.5 


J1909+0749 


4.39 


12.8 


35.6 


J1916+1225 


5.19 


12.4 


34.9 


J1917+1737 


7.21 


11.5 


32.5 


J1919+1314 


6.38 


12.2 


32.9 


J1919+1745 


7.29 


12.3 


30.9 


J1924+1631 


5.11 


13.5 


32.7 


J1928+1746 


4.92 


12.0 


36.2 


J1934+2352 


4.34 


12.7 


36.0 


J1938+2010 


6.51 


12.2 


32.6 


J1940+2337 


5.05 


12.8 


34.3 


J1941+2525 


5.36 


13.3 


32.7 


J1946+2535 


6.16 


12.2 


33.2 


J1947+1957 


6.68 


11.5 


33.7 


J1948+2551 


5.54 


12.1 


34.7 


J1949+2306 


8.23 


11.6 


30.3 


J1953+2732 


7.08 


12.2 


31.5 


J2005+3547 


7.54 


11.6 


31.7 


J2006+3102 


5.02 


12.3 


35.3 


J2007+3120 


5.79 


12.5 


33.4 


J2009+3326 


7.19 


12.2 


31.3 


J2010+2845 


7.99 


11.4 


31.3 


J2010+3230 


6.80 


12.4 


31.7 


J2011+3331 


6.92 


12.1 


31.9 


J2013+3058 


7.46 


11.3 


32.4 


J2018+3431 


6.52 


11.9 


33.1 



10 



D. J. Nice et al. 



Table 3 

Galactic Positions, Distance, and Scattering Properties. 

NE2001 Predictions'' 



Measured Values'' 



PSR 


I 


b 


DM 


d 


TlGHz 


TlGHz 


ADM 


'^comp 


xi 




(°) 


(°) 


(pccm~^) 


(kpc) 


(ms) 


(ms) 


(pccm~^) 






J1924+1631 


51.405 


0.318 


518.5 


14.0 


0.50 


49(2) 


-10.4(8) 


2 


1.3 


J1858+0346 


37.084 


0.182 


386.6 


7.0 


4.97 


42(3) 


-3.0(4) 


3 


1.1 


J2010+3230 


70.391 


-0.495 


371.8 


12.0 


0.08 


24.8(6) 


-3.8(1) 


1 


1.0 


J1953+2732 


64.205 


0.059 


195.4 


7.0 


0.01 


21(3) 


-2.8(2) 


3 


1.1 


J1919+1314 


47.895 


-0.088 


613.5 


15.2 


1.87 


18.7(9) 


-2.11(5) 


2 


1.5 


J1856+0245 


36.008 


0.057 


623.5 


9.0 


23.31 


18.0(6) 


-3.0(1) 


2 


1.0 


J1855+0205 


35.281 


0.007 


867.3 


11.6 


57.33 


16.3(9) 


-2.1(1) 


2 


1.4 


J2005+3547 


72.581 


2.180 


401.7 


>50.0 


0.07 


11(1) 


-1.85(9) 


2 


1.5 


J1850+0423 


36.718 


2.231 


265.9 


6.4 


0.28 


7.8(8) 


0.0(2) 


2 


1.3 


J1941+2525 


61.037 


1.263 


314.4 


9.8 


0.03 


7.4(6) 


-2.4(2) 


2 


1.5 


J1948+2551 


62.207 


0.131 


289.3 


8.8 


0.02 


6.3(2) 


-0.83(1) 


4 


1.5 


J1949+2306 


59.931 


-1.423 


196.4 


7.0 


0.01 


5.4(7) 


-1.0(1) 


3 


1.4 


J1916+1225 


46.811 


0.225 


265.3 


7.1 


0.59 


5.3(2) 


-0.155(9) 


4 


1.6 


J1946+2535 


61.809 


0.283 


248.8 


8.1 


0.01 


4.54(7) 


-0.97(2) 


3 


1.8 


J1938+2010 


56.115 


-0.672 


327.7 


9.5 


0.05 


4(1) 


0.1(4) 


2 


1.1 


J2011+3331 


71.320 


-0.049 


298.6 


8.8 


4.16 


4.39(9) 


-0.84(2) 


3 


1.9 


J1940+2337 


59.397 


0.529 


252.2 


8.2 


0.02 


4.2(5) 


-0.8(1) 


2 


1.3 


J2007+3120 


69.042 


-0.535 


191.6 


6.8 


0.01 


3.8(8) 


-0.52(6) 


2 


1.6 


J1919+1745 


51.897 


1.987 


142.3 


5.2 


0.01 


2.9(2) 


-0.60(5) 


4 


1.8 


J1909+0749 


41.909 


-0.340 


539.4 


9.5 


6.31 


1.9(1) 


-0.32(4) 


5 


1.5 


J1904+0738 


41.180 


0.680 


278.3 


6.3 


1.01 


1.61(7) 


-0.30(2) 


3 


1.4 


J2018+3431 


73.044 


-0.841 


222.3 


7.2 


0.02 


1.2(1) 


-0.35(3) 


2 


1.3 


J1905+0902 


42.555 


1.056 


433.5 


9.1 


1.78 


0.9(2) 


-0.22(6) 


3 


1.6 


J2013+3058 


69.485 


-1.886 


148.8 


6.0 


0.00 


<0.9 


-0.10(6) 


3 


1.6 


J1917+1737 


51.528 


2.418 


208.1 


7.0 


0.02 


<1.0 


-0.1(3) 


2 


1.2 


J1928+1746 


52.931 


0.114 


176.7 


5.8 


0.04 


<0.6 


-0.11(4) 


5 


1.3 


J1848+0351 


36.058 


2.366 


336.7 


7.9 


0.41 


<1.0 


0(1) 


1 


1.5 


J2006+3102 


68.667 


-0.530 


107.2 


4.7 


0.00 


<1.0 


0.03(6) 


5 


1.4 


J2009+3326 


71.103 


0.124 


263.6 


7.9 


1.73 


<1.0 


0.1(1) 


3 


1.3 


J1934+2352 


58.966 


1.814 


355.5 


11.6 


0.03 


<1.0 


0.1(2) 


2 


1.6 


J2010+2845 


67.210 


-2.472 


112.5 


5.0 


0.00 


<1.0 


-0.07(3) 


3 


3.5 


J1909+0641 


40.941 


-0.940 


36.7 


2.2 


0.00 


<1.0 


0.00(3) 


3 


1.6 


J0540+3207 


176.719 


0.761 


62.0 


1.8 


0.00 


<1.0 


-0.060(5) 


2 


2.6 


J0628+0909 


202.190 


-0.851 


88.3 


2.5 


0.00 


<1.0 


0.144(9) 


1 


3.7 


J1947+1957 


56.987 


-2.658 


185.9 


6.8 


0.01 


<1.0 


-0.03(4) 


2 


1.5 



Values predicted based on /, b. and DM, usin g the NE2001 electron density model of |Cordes &€ Lazio| ( |2002[ |. 
Values measured from fits to pulsar profiles (^6.2[l. Figures in parentheses arc uncertainties in the last digit quoted. 
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Table 4 

Flux Densities, Spectral Indices, and Luminosities. 



PSR 



•S'l400 
(mjy) 



(Jykpc2) 



J0540+3207 
J0628+0909 
J1848+0351 
J1850+0423 
J1855+0205 
J1856+0245 
J1858+0346 
J1904+0738 
J1905+0902 
J1909+0641 
J1909+0749 
J1916+1225 
J1917+1737 
J1919+1745 
J1919+1314 

J1924+1631 
J1928+1746 
J1934+2352 
J1938+2010 
J1940+2337 
J1941+2525 
J1946+2535 
J1947+1957 
J1948+2551 
J1949+2306 
J1953+2732 
J2005+3547 
J2006+3102 
J2007+3120 
J2009+3326 
J2010+2845 
J2010+3230 
J2011+3331 
J2013+3058 
J2018+3431 



0.34(2) 


-2.0(3) 


0.058(3) 


-0.7(5) 


0.066(2) 


-2.1(4) 


0.199(4) 


-2.9(2) 


0.193(5) 


-1.4(3) 


0.58(2) 


-0.4(2) 


0.190(9) 


-3.0(4) 


0.23(1) 


-1.1(4) 


0.097(6) 


-1.7(5) 


0.111(3) 


-1.9(2) 


0.226(8) 


-1.3(3) 


0.094(2) 


-0.6(2) 


0.047(1) 


-1.4(4) 


0.190(8) 


-2.6(4) 


0.224(6) 


-1.3(3) 


0.088(4) 


-0.9(3) 


0.278(8) 


0.0(3) 


0.062(2) 


-3.4(3) 


0.103(9) 


-1.9(7) 


0.069(4) 


0.1(7) 


0.240(8) 


-0.8(3) 


0.480(7) 


-1.6(2) 


0.081(6) 


-1.5(8) 


0.622(7) 


-0.7(1) 


0.097(4) 


-3.3(3) 


0.106(7) 


-3.0(6) 


0.252(9) 


-2.2(4) 


0.27(1) 


-1.8(4) 


0.123(5) 


-2.9(4) 


0.15(1) 


-4.1(6) 


0.43(2) 


-2.1(3) 


0.118(5) 


-3.5(3) 


0.384(8) 


-2.5(2) 


0.067(5) 


-3.0(7) 


0.247(9) 


-1.2(3) 



0.384 


1.1 


0.571 


0.4 


0.665 


4.1 


0.687 


8.1 


1.313 


24.5 


0.602 


45.5 


0.494 


10.2 


0.029 


9.1 


0.263 


8.5 


1.362 


0.5 


0.098 


1.0 


1.863 


4.4 


0.004 


2.3 


0.174 


5.5 


2.378 


53.9 


3.794 


17.6 


1.866 


9.5 


1.185 


8.4 


0.233 


8.8 


1.911 


5.0 


0.097 


22.6 


4.376 


30.9 


0.027 


3.4 


5.619 


47.9 


1.089 


7.8 


0.748 


9.5 


0.082 


a 


0.157 


6.4 


0.291 


6.1 


0.261 


19.7 


0.373 


9.9 


0.902 


31.6 


1.215 


32.9 


0.253 


2.4 


0.693 


13.1 



Unknown L1400 due to unknown distance. 
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Table 5 


















Pulse Widths 










PSR'' 


P (ms) 


1170 


Wm (ms) 
1370 1470 


1570 




W,, 


o/P 






1170 


1370 


1470 


1570 






MHz 


MHz 


MHz 


MHz 


MHz 


MHz 


MHz 


MHz 


0540+3207 


524.3 


11.0 


11.7 


12.0 


12.4 


0.021 


0.022 


0.023 


0.024 


0628+0909 


1241.4 


7.6 


7.7 


8.0 


8.4 


0.006 


0.006 


0.006 


0.007 


1848+0351 


191.4 


9.9 


9.9 


10.3 


12.0 


0.052 


0.052 


0.054 


0.063 


1850+0423 


290.7 


24.6 


24.3 


24.7 


24.6 


0.084 


0.084 


0.085 


0.085 


1855+0205 


246.8 


20.5 


9.4 


7.6 


6.6 


0.083 


0.038 


0.031 


0.027 


1856+0245 


80.9 


23.8 


17.6 


14.6 


14.2 


0.295 


0.218 


0.180 


0.175 


1858+0346 


256.8 


41.4 


33.9 


36.7 


32.3 


0.161 


0.132 


0.143 


0.126 


1904+0738 


209.0 


3.2 


3.0 


2.9 


2.8 


0.015 


0.014 


0.014 


0.013 


1905+0902 


218.3 


3.9 


3.7 


3.9 


3.9 


0.018 


0.017 


0.018 


0.018 


1909+0641 


741.8 


8.2 


7.9 


7.9 


7.6 


0.011 


0.011 


0.011 


0.010 


1909+0749M 


237.2 


5.6 


4.7 


5.2 


5.0 


0.024 


0.020 


0.022 


0.021 


1909+07491 


237.2 


4.0 


3.7 


3.7 


3.4 


0.017 


0.016 


0.016 


0.014 


1916+1225 


227.4 


2.0 


1.7 


1.7 


1.6 


0.009 


0.008 


0.007 


0.007 


1917+1737 


334.7 


6.1 


5.5 


5.9 


5.7 


0.018 


0.016 


0.018 


0.017 


1919+1314 


571.4 


15.6 


13.2 


11.5 


11.6 


0.027 


0.023 


0.020 


0.020 


1919+1745 


2081.3 


63.6 


62.4 


62.4 


62.5 


0.031 


0.030 


0.030 


0.030 


1924+1631 


2935.2 


35.1 


27.1 


27.2 


24.1 


0.012 


0.009 


0.009 


0.008 


1928+1746 


68.7 


3.4 


3.5 


3.4 


3.4 


0.050 


0.051 


0.050 


0.049 


1934+2352 


178.4 


4.5 


3.7 


4.2 


3.3 


0.025 


0.021 


0.024 


0.018 


1938+2010 


687.1 


24.7 


24.5 


25.0 


24.7 


0.036 


0.036 


0.036 


0.036 


1940+2337 


546.8 


10.6 


10.4 


9.4 


10.3 


0.019 


0.019 


0.017 


0.019 


1941+2525 


2306.2 


33.5 


36.6 


35.6 


33.4 


0.015 


0.016 


0.015 


0.015 


1946+2535 


515.2 


9.5 


8.7 


8.6 


8.5 


0.018 


0.017 


0.017 


0.016 


1947+1957 


157.5 


6.9 


6.1 


6.5 


5.9 


0.044 


0.039 


0.041 


0.038 


1948+2551 


196.6 


9.1 


9.0 


9.2 


9.0 


0.046 


0.046 


0.047 


0.046 


1949+2306 


1319.4 


17.8 


17.8 


19.6 


21.5 


0.013 


0.014 


0.015 


0.016 


1953+2732 


1334.0 


29.7 


29.0 


26.3 


26.7 


0.022 


0.022 


0.020 


0.020 


2005+3547 


615.0 


23.3 


19.6 


20.8 


18.9 


0.038 


0.032 


0.034 


0.031 


2006+3102 


163.7 


12.2 


11.8 


11.8 


11.8 


0.074 


0.072 


0.072 


0.072 


2007+3120 


608.2 


8.7 


8.9 


8.7 


8.8 


0.014 


0.015 


0.014 


0.015 


2009+3326 


1438.4 


24.5 


25.8 


26.7 


29.7 


0.017 


0.018 


0.019 


0.021 


2010+2845 


565.4 


14.6 


14.1 


13.6 


13.5 


0.026 


0.025 


0.024 


0.024 


2010+3230 


1442.4 


30.6 


28.1 


28.8 


26.8 


0.021 


0.019 


0.020 


0.019 


2011+3331 


931.7 


30.7 


30.6 


30.6 


30.1 


0.033 


0.033 


0.033 


0.032 


2013+3058 


276.0 


3.0 


2.8 


2.9 


3.0 


0.011 


0.010 


0.010 


0.011 


2018+3431 


387.7 


6.8 


6.6 


6.7 


6.5 


0.017 


0.017 


0.017 


0.017 



Main pulse and intcrpulsc denoted by M and I, respectively. 



Timing and Scattering of PALFA Pulsars 



13 



•Ph - 

GO 
O 
J - 





:i 


MM 


1 1 i^yii 


1 1 1 M MM 1 1 III 1 M L III 


III 


U- ' J 


11 


=— 




c 


'^^><^^^ '• "^^<^ 


• 




12 


__ 


\.^^ 


^^ 






_: 


13 
14 


''^.,,^^ 


vo^^' 






• 
• 

•• 
• • 

• 


• 




I 


a* ^^ 


^0/7 


^>^gM|H||^HBMBi^; * 




^^^^ I 


15 


= ~, 


vo 


c 


><r^^^^^^^^\ 


• 


~= 


16 


E- 


^^\ 


^^^/-^" 


^""s^^rf '^* *. * 


^--^ 


\J 


17 


=- 


vo^^' 






^ 


-= 


18 


E- 


"^"^^--^ 


^^^-^^'^ 




^^ 


\^-l 


19 


1^ 


vo^ t 








-= 




=- 


• 
• 


• •• 

• • 


• ^""■---.,.,.^^ 


^^^ 


^-- 


"^^^^^ 


20 

21 


^^ 


O This paper 
A Other PALFA 






~ 








• ATNF catalog 




~ 


22 


















(>jf>j 


= 1 


1 1 III 


1 


1 1 1 III 


1 II 


">-d 



0.001 



0.01 



0.1 
P (s) 



10 



Figure 1. Pulsar rotation periods and period derivatives. Circles (blue in online version) indicate pulsars in the present paper; triangles 
(blue in online version) indicate other published pulsars discovered in the PALFA survey; and dots (red in online version) are all known 
pulsars in the Galactic disk (Hobbs fc Manchester, 2013;) . 
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Figure 2. Residual pulse arrival times of all pulsars in this paper are given in Figures[2]through|6] The points are daily average residuals 
at Arecibo (circles, blue in online version) and JodrcU Bank (x's, red in online version). For each pulsar, the residual plot was made 
by performing a timing fit, possibly including wrhitening and glitch parameters as described in the text, to find the best-fit parameters; 
freezing the position at its best-fit value; removing whitening and glitches from the model; re-fitting for spin period and spin-down rate; 
and averaging the resulting residual pulse arrival times within each day at each observatory. 
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Figure 3. See the caption for Figure [2] 
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Figure 4. See the caption for Figure [2] 
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Figure 5. See the caption for Figure [2] 
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Figure 6. See the caption for Figure [2] 
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Figure 7. Pulsar positions projected onto the Galactic plane. Pulsar positions are based on the coordinates and distances given in Tablels] 
The position of the Sun is indicated by a small circle and dot at (x, y) = (0, 8.5 kpc). Straight lines emerging from the Sun indicate regions 
of low Galactic latitude within t he declination ra nge of the Arecibo telescope. Curved lines indicate Galactic spiral arms according to the 
NE2001 electron density model {Cordes & Lazio 2002|)). Circles (blue in online version) indicate pulsars in the present paper; triangles 
(blue in online version) indicate other published pulsars discovered in the PALFA survey; and dots (red in online version) are all known 
pulsars in the Galactic disk IjHobbs fc Manchester|2013^. 
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Figure 8. Profiles of the pulsars in each of the four observational subbands are given in Figures [s] and |9] Each profile is a full pulse 
period. Profiles are sums of all Arccibo data records which produced valid TOAs for a given pulsar at a given frequency. Da ta r eco rds 
were aligned using the final timing solution for each pulsar. Plots are total intensity (summed polarizations). See Figures |10| and |ll| for 
expanded views of the pulse peaks. 
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Figure 9. See caption to figure [8] 
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Figure 10. Closeup of the pulse peak of each pulse profile in Figure Is] Horizontal range varies from pulsar to pulsar depending on the 
pulse width and is given above each set of profiles, e.g., "0540+3207(10%)" means that the horizontal extent of the profiles of 0540+3207 
is 10% of the pulse period. Two plots are shown for 1909+0749, labeled by "M" and "I," for main pulse and interpulse; these correspond 
to the leading ("M") and trailing components ("I") of the profile shown in Figure [s] These two components are similar in total energy, and 
the designation of "M" and "I" is arbitrary. They are plotted on the same verticaTscale, and the centers of the two plots are separated by 
half a pulse period. 
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Figure 11. Closcup of the pulse peak of each pulse profile in Figurcl9] See eaption to Figure [To] for further details. 
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Figure 12. Sample pulsar scattering model fits. Top: PSR J1848+0351, PSR J1855+0205; Middle: PSR J1938+2010, PSR J1948+2551; 
Bottom: PSR J1953+2732, PSR J2005+3557. For each pulsar, closeups of the pulse peak are shown for the four observing subbands 
(top to bottom: 1570, 1470, 1370, and 1170 MHz). Superimposed on each profile are the profile models for that frequency both with 
scattering (black line) and without scattering (gray line). Of the pulsars shown here, J1855+0205 and J1953+2732 show strong scattering; 
J1848+0351 and J1938+2010 show almost no scattering; and J1948+2551 and J2005+3547 have high covariance between the Gaussian 
profile model and the scattering tail. 
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Figur e 13. Scattering time scales of tlie pulsars, comparing our measured time scales with the predictions of the NE2001 electron density 
model | |Cordes fc Lazio|2002| |. 
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Figure 14. Left. Solid contours (green in online version): HESS 18 57+026 TeV source. Dashed contours (magenta in online version): 
Fermi Large Area Telescope source coincident with HESS J1857+026 ( [Rousseau et al.'2012'). Contours correspond to test statistic source 
significance levels of 25, 16, and 9 (5, 4, and 3 a). Dashed rectangle (red in online version): XMM-Newton EPIC field of view. Dashed 
square (blue in online version): Chandra ACIS-I 1-8 keV image. Black square: region shown in right image; pulsar position is in the center. 
Right. Chandra ACIS-I image in the 1-8 keV band centered on PSR J1856-I-0245. The ellipse (red in online version) shows the position 
uncertainty from the pulsar timing solution (Table [TJ. Pixel randomization was removed from the pipeline processing. 



